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ABSTRACT 

The X-ray properties of twenty ~lMyr old O, B, and A stars of the Orion 
Trapezium are examined with data from the Chandra Orion Ultradeep Project 
(COUP). On the basis of simple theories for X-ray emission, we define two classes 
separated at spectral type B4: hotter stars have strong winds that may give rise to 
X-ray emission in small- or large-scale wind shocks, and cooler stars that should 
be X-ray dark due to their weaker winds and absence of outer convection zones 
where dynamos can generate magnetic fields. Emission by late-type magnetically 
active companions may be present in either class. 

Sixteen of the 20 stars are detected with a wide range of X-ray luminosi- 
ties, logL x [ergs -1 ] ~ 29 — 33 and X-ray efficiencies log(L x /L bol ) ~ —4 to —8. 
Only two stars, 8 1 OriD (BO. 5) and NUOri (Bl), show exclusively the constant 
soft-spectrum emission at log (L x /Lt, i) ~ —7 expected from the standard model 
involving many small shocks in an unmagnetized radiatively accelerated wind. 
Most of the other massive 07-B3 stars exhibit some combination of soft-spectrum 
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wind emission, hard-spectrum flaring, and/or rotational modulation indicating 
large-scale inhomogeneity. Magnetic confinement of winds with large-scale shocks 
can be invoked to explain these phenomena. This is supported in some cases by 
non-thermal radio emission and/or chemical peculiarities, or direct detection of 
the magnetic field (^ 1 OriC). 

Most of the stars in the weak-wind class exhibit X-ray flares and logL x < 
31 erg s _1 , consistent with magnetic activity from known or unseen low-mass 
companions. In most cases, the X-ray spectra can be interpreted in terms of a 
two-temperature plasma model with a soft component of 3 — 10 MK and a hard 
component up to 40 MK. All non-detections belong to the weak-wind class. A 
group of stars exhibit hybrid properties - flare-like behavior superimposed on 
a constant component with logL x ~ 32 erg s _1 - which suggest both magnetic 
activity and wind emission. 

Subject headings: open clusters and associations: individual: (Orion nebula clus- 
ter) - binaries: general - stars: early-type - stars: pre-main sequence - stars: 
winds, outflows - X-rays: stars 

1. Introduction 

Among normal stars, X-ray emission is often strongest in the youngest stellar popula- 
tions. X-ray emission from OB stars was first reported from observations with the Einstein 
observatory (Harnden et al. 1979; Pallavicini et al. 1981), and dozens were detected with the 
ROSAT All-Sky Survey with characteristic efficiency log (L x /L bol ) ~ —7 (Berghofer et al. 
1997). Their soft and nearly constant emission was attributed to the integrated emission 
from a myriad of small shocks due to instabilities in the radiatively-driven outflow (see Lucy 
& White 1980; Owocki & Cohen 1999; Kudritzki & Puis 2000, for a review of the wind 
properties of massive stars). Their soft spectrum and mildly variable X-ray emission stood 
in contrast to the hard and often wildly variable X-ray emission of lower mass T Tauri stars. 
The variability of these latter stars is attributed to flares resulting from reconnection events 
in large magnetic loops as in the Sun and other magnetically active stars (Feigelson & Mont- 
merle 1999; Favata & Micela 2003). Thus, early- type stars were viewed to radiate X-rays via 
a thermal plasma produced by a thermal heating mechanism (wind shocks), while late-type 
stars radiate via a thermal plasma produced by a non-thermal heating mechanism (magnetic 
reconnection flares). 

Recent studies by telescopes with sensitivity to hard X-rays and high-resolution spectral 
capabilities have thrown this simple view into confusion. The brightest star in the Orion 
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Trapezium cluster, the 07 star 9 1 Ori C, exhibits a hard spectral component and a strong 
rotational modulation inconsistent with spherically symmetric wind emission. The Trapez- 
ium 09.5 star 6 2 Ori A showed a powerful flare superposed on wind-like emission (Feigelson 
et al. 2002). Spectral lines in these and other O stars are broadened indicating production 
in the wind acceleration zone, but the line strengths, shapes and absorption often differ from 
those predicted by the standard wind-shock model (e.g. Miller et al. 2002; Schulz et al. 2003; 
Cohen et al. 2003). Some of these effects can be attributed to the confinement, and result- 
ing large-scale shocks, of the wind by a strong stellar magnetic field (Babel & Montmerle 
1997a; Ud'Doula & Owocki 2002; Schulz et al. 2003). This is supported by the detections of 
nonthermal radio emission and atmospheric chemical peculiarities in some of these stars. 

The X-ray properties of intermediate-mass late B and A stars have also been confusing, 
but for different reasons. With neither powerful fast winds, nor outer convection zones, no 
X-rays were expected. Nonetheless, a substantial number are detected with a wide range of 
X-ray luminosities (e.g. Schmitt et al. 1985; Caillault & Zoonematkermani 1989; Simon et al. 
1995; Cohen et al. 1997). In some cases, this emission can be attributed to the presence 
of low-mass, late-type star companions to the intermediate-mass stars. This was directly 
demonstrated by resolving several visual binaries with the Chandra X — ray Observatory 
(Stelzer et al. 2003; Behar et al. 2004). But other cases are less clear: some X-ray sources 
appear to coincide precisely with the B star (Stelzer et al. 2005), and the youngest A and 
B stars (Herbig Ae/Be stars) sometimes outshine their lower mass companions (Feigelson 
et al. 2003), perhaps due to star-disk magnetic interactions (Hamaguchi et al. 2005). 

These issues can be fruitfully addressed in the Orion Nebula Cluster (ONC) where a 
roughly coeval population of young stars, covering a variety of early spectral types, can be 
observed in a single X-ray exposure. With its high spatial resolution and excellent astrometry, 
Chandra is the best instrument for a study of this rich and crowded field. Results from 
Chandra observations carried out in 1999 and 2000 already showed that most OBA stars 
in the ONC are X-ray emitters 1 above logLf ~ 28 ergs -1 (Feigelson et al. 2002; Flaccomio 
et al. 2003). The luminosities and variability of some OB stars are consistent with wind 
emission. Similar to findings from earlier X-ray observations, the emission levels are seen 
to scatter around the canonical L x /Lb i — 10~ 7 relation. This relation has no immediate 
physical interpretation but could be recovered in a scaling analysis taking account of the 
detailed wind-shock structure and wind opacity (Owocki & Cohen 1999). The large scatter 
of the data around the canonical value has remained unexplained, but may be related to 



following Getman et al. (2005a) and other COUP studies, L t is the X-ray luminosity in the 'total' 0.5 — 8 
keV band assuming a distance of 450 pc to the ONC. We later refer to L s in the 'soft' 0.5 — 2 keV band, 
and L t . c where the total band luminosity is corrected for absorption. 
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the fact that L x and Lboi are coupled to each other only indirectly via the wind density 
parameter Mq/v^. Contrary to expectation, the ~ 30 M massive Trapezium star 8 2 Ori A 
exhibited rapid high-amplitude variability characteristic of magnetic flares (Feigelson et al. 
2002). Together with other works on early- type stars, these results hint that the X-ray 
emission mechanism in hot stars is more complex than the simple wind-shock picture outlined 
above, suggesting in particular that transient magnetic phenomena may be important even 
in massive stars. In this paper, we seek to clarify these issues and elucidate the respective 
roles of wind shocks and magnetic activity in the X-ray emission of early-type stars. 

Our study of early-type stars in the ONC is based on the nearly-continuous 9.7 day 
observation of the Orion Nebula made with Chandra in 2003. Known as the Chandra Orion 
Ultradeep Project (COUP), this is the most comprehensive X-ray study of a young stellar 
cluster performed to date. The observation, data analysis, source lists and characteristics 
are described in detail by Getman et al. (2005a). After describing the sample and presenting 
its basic X-ray properties (§2), we examine their X-ray spectra (§3) and various types of 
variability (§4). Section 5 with Appendix A is dedicated to multiple systems, §6 to systems 
with disks, and Appendix B to a new X-ray detection of an early-type Trapezium star in a 
crowded region of the field. In the discussion (§7), we confront different possible emission 
mechanisms for our targets. 



2. Early-type stars in the COUP 

2.1. Classification 

In order to investigate whether or not magnetic phenomena are involved in the X-ray 
emission mechanisms for the early-type stars of the ONC, we first divide the whole COUP 
sample of ONC members confirmed by optical spectroscopy ('optical' COUP sample) into 
three groups distinguished by the simplest emission mechanism associated as a class with 
their spectral types (§1): the hottest stars dominated by wind shocks, the cooler stars 
dominated by magnetic activity, and an intermediate 'blank zone'. In this last class no X- 
rays are expected, because the stars have no, or only a very shallow (outer) convective zone 
(at most a few percent of the stellar radius; Siess et al. 2000), hence no magnetic activity, 
and have weak, slow radiative winds or no wind at all. We seek a quantitative basis for 
dividing these classes here. 

In hot massive stars, a radiative force is exerted on heavy atoms in their atmospheres by 
photospheric ultraviolet photons, transmitting outward momentum to the entire atmosphere 
via Coulomb collisions (e.g., Kudritzki & Puis 2000). When the radiative acceleration g ra d 
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of the gas exceeds the local gravity g*, a wind is triggered with net acceleration g ra d — g±. 
The wind velocity increases up to the sonic point, reaching a ballistic regime with "terminal 
velocity" > v csc , where v esc is the escape velocity of the star. The development of a 
radiatively accelerated wind thus depends both on a star's radiation field indicated by its 
spectral type (or T c g) and its gravity g*. Metallicity should also play a role but is not 
considered here, since all stars in our sample are assumed to have the same metallicity. 

Early models for radiatively driven winds (e.g., Castor et al. 1975; Abbott et al. 1979) 
gave a simple, continuous relationship between v^, the mass-loss rate M, and T e g, so that 
the transition between slow winds and fast winds was gradual. However, observations (par- 
ticularly in the X-ray band) suggest a sharper transition both in Voo and M around mid-B 
spectral types. Taking into account the Doppler shifts and shadowing of the photospheric 
lines, Babel (1996) calculated the locus of the stationary wind regime in the two-dimensional 
g — T cS plane. He found that for low temperatures (e.g., T eff < 10, 000 K), only weak-gravity 
stars can have winds, whereas on the main-sequence, for high temperatures (T e ff > 20, 000 K), 
winds can exist even for high-gravity stars. In practice, the transition between the weak-wind 
and fast- wind regimes near the main-sequence takes place for T e s ~ 14, 000 — 18, 000 K; i.e., 
between spectral types B3 and B5 on the temperature scale of Hillenbrand & White (2004). 

This choice of boundary between strong and weak winds is supported by two lines of 
evidence. First, the dependence of Voo/v csc on T c fj saturates at foo/^csc ~ 3 for T e g > 20, 000 K 
(Kudritzki & Puis 2000). Second, a simple model of X-ray emission from wind shocks shows 
a break from L x oc L\ ol to L x oc L\ ol around mid-B stars (Owocki & Cohen 1999), basically 
due to a rapid change in the opacity of the winds to UV and X-ray photons. 

We therefore classify our sample in two groups based on their spectral type range and 
expected source of X-ray emission: 

1. stars with spectral type O through B3 that have strong, optically thick winds (SW), 

2. stars with spectral type B5 through A9 that have weak, optically thin winds (WW), 

For comparison, we will also make use of the properties of cool, convective stars with spectral 
type F0 and later that are magnetically active (MA). This physical classification will help 
us identify and interpret deviations from the simple model predictions outlined in §1. 

The coincidence of the lower boundary of the WW stars with the appearance of outer 
convection zones in cool stars is demonstrated in Fig. 1. It shows the Hertzsprung- Russell 
diagram with evolutionary interior models by Siess et al. (2000). The jagged curves denote 
the extent of the convective envelope as a fraction of the stellar radius. The data points 
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show the 'optical COUP sample' of Hillenbrand (1997) (henceforth H97). In all WW stars 
(shown as large circles), the convective envelopes comprises at most 5 % of the stellar radius. 

Roughly two thirds (21/32) of the OBA stars in the ONC are in the field-of-view of the 
COUP. One star, JW602, has poorly characterized optical properties (the spectral type is 
classified as 'B:') and is discarded. For the remaining 20 stars, 9 fall in the SW sample and 
11 in the WW sample 2 . 

Table 1 summarizes stellar parameters for the 20 COUP OBA stars, in order of decreas- 
ing effective temperature. We adopt the V magnitude, Ay, and L bol from Getman et al. 
(2005a), who had extracted the parameters from the original data by H97. When a range 
of spectral types are present in the literature, a spectral type was assigned corresponding 
to an estimated average of the group (Hillenbrand, unpublished) and converted to an effec- 
tive temperature using the scale of Hillenbrand & White (2004). Given the possibility of 
yet unresolved MA late-type companions, we searched the literature for multiplicity of our 
targets which is summarized in columns 7 — 9. Finally, in order to investigate the possibility 
of magnetic star-disk interactions, as in young protostars and T Tauri stars (e.g., Montmerle 
et al. 2000; Bouvier et al. 1999), we provide in column 10 flags for stars with infrared (IR) 
excess indicating the presence of circumstellar material. 



2.2. X-ray properties 

The time- averaged X-ray properties for the 20 stars integrated over the 9.7-day COUP 
exposure are presented in Table 2. Except for the hardness ratio (column 9), the values are 
reproduced from the tables of Getman et al. (2005a) who give details of their derivation. The 
first three columns give the COUP sequence number, the optical counterpart from Table 1, 
and the spectral type. Following columns provide the off-axis angle and the offset between 
X-ray source and the near-IR stellar positions. The net source counts in the broad energy 
band (0.5 — 8.0 keV) after background subtraction and upper limits for undetected sources 
are given in column 6. 

Four of the twenty targets are heavily piled-up in the ACIS CCD detector: COUP 809 
(^OriC), 1232 (# 2 OriA), 745 (^OriA), and 732 (^OriE). In these cases, only 1 -4% 
of the photons incident on the detector are extracted from an annular region around the 



2 Two of the WW stars with late-B spectral type have unknown membership probabilities: 8 1 OriE and 
9 1 Ori A. Since B stars are intrinsically rare, we consider it unlikely that these stars are not members of the 
ONC. 
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piled-up core of the point spread function (PSF). The spectral and luminosity analysis takes 
this into account in a self-consistent manner. The luminosity values obtained for heavily 
piled-up sources are subject to systematic uncertainties of order ±20 %. The treatment of 
piled-up COUP sources is described in §6 of Getman et al. (2005a). 

All O to mid-B (SW) stars are detected, including the B2 star Par 1772, a faint X- 
ray source in a crowded region which had not been detected in previous Orion Nebula 
X-ray studies (see Appendix A and Fig. 10 for details). Four mid-B to A (WW) stars are 
undetected: ^OriF (B8), JW531 (Al), JW108 (A3), and JW608 (A5). None of these 
stars has exhibited X-ray emission in previous observations. Upper limits to their X-ray 
luminosities were computed with XSPEC from the counts given in Table 2, assuming an 
isothermal plasma with energy kT = 1 keV, abundance of 0.3 solar, and absorbing column 
derived from the visual absorptions given in Table 1 according to N H [cm~ 2 ] = 1.6 10 21 
A y [mag] (Vuong et al. 2003). 

We also list the PSF fraction included in the photon extraction area (column 7), and the 
effective exposure time (column 8) at that location in the ACIS detector. The PSF fraction 
is ~ 90 % for most COUP sources, but is lower when the source suffers from pile-up or has 
a very close neighbor. In such cases, special photon extraction regions have been defined. 

The hardness ratio in column 10 is defined as H/S, where H is the number of counts in 
the hard band (E > 1.8 keV) and S is the number of counts in the soft band (E < 1.8 keV). 
This differs from the hardness ratios defined by Getman et al. (2005a). Our choice of the 
boundary between the hard and soft band is near to the peak of the observed spectrum of our 
targets, and therefore guarantees a significant fraction of the total collected source photons 
in each of the bands. In the last column, we provide the number of the respective X-ray 
source in the previous Chandra ACIS observation (Feigelson et al. 2002, hereafter FBG02). 

3. Spectral Analysis 

We adopt the procedure to extract the X-ray photons described by Getman et al. 
(2005a), but we conduct a more detailed spectral analysis than their automated fitting 
procedure, which did not take in account visual absorptions, unusual abundance or temper- 
ature distributions, or variability. In addition, as illustrated in their Figure 6, their fits may 
give incorrect absorption column densities. Our resulting spectral fits to thermal plasma 
models are summarized in Table 3. X-ray fluxes and luminosities were computed from the 
spectral fits using the XSPEC flux command assuming a distance of 450 pc. We list both 
the absorbed (L t ) and the absorption-corrected (£ t ,c) luminosities in the 0.5 — 8 keV broad 
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band. 

Spectral modeling was first performed based on a one- or two-temperature plasma in 
collisional ionization equilibrium using the MEKAL code (Mewe et al. 1985, 1995). Elemental 
abundances were fixed at Z — 0.3 Z with abundance ratios matching those of Anders & 
Grevesse (1989). This choice of subsolar abundances follows Getman et al. (2005a), and 
has been proven appropriate in X-ray studies of various star forming regions (e.g. Nakajima 
et al. 2003; Feigelson et al. 2002; Preibisch 2003). The spectra of over half of the ONC early- 
type stars are well-fit by an absorbed, two-temperature (2-T) plasma subject to soft X-ray 
absorption by a cold intervening medium with column density Ah- As noted by Getman 
et al., in some sources the spectral models are not unambiguous: a large uncertainty is 
introduced by the choice of A H , resulting in either very cool and heavily absorbed plasma 
or somewhat hotter plasma subject to light absorption. We have also tested models with 
A H fixed on the value expected from the Vuong et al. (2003) extinction law, but found after 
visual inspection of the best fits that the models with free column density are more adequate. 

In about half of the sample, the spectral parameters resulting from our analysis are 
almost identical to those presented in Getman et al. (2005a). Discrepancies for the remaining 
cases can be attributed to two causes. First, Getman et al. sometimes give parameters 
of a 1-T model while, to achieve a uniform description of our sample, we prefer the 2-T 
models. The exception is COUP 995 where the 2-T fits result in lower temperatures and 
unreasonably high luminosity. Second, the fitting procedure encounters a local rather than 
global minimum in one of the studies. In these cases, we retain the model which both gives 
a satisfactory fit (Xrcd ~ 1) anc ^ nas absorption consistent to expected ONC values. That 
is, the derived column density is neither extremely low (log Ah [cm' 2 ] < 20) nor extremely 
high (log A H [cm- 2 ] > 22). 

In 4 of the brightest COUP hot stars, no 2-T model provides an adequate description of 
the observed spectrum. COUP 732 (^OriE), 745 (^OriA), and 809 (^OriC), which all 
suffer strongly from photon pile-up in the ACIS detector, show a high-energy excess which 
calls for the inclusion of a third high-temperature component in the model. This high- 
temperature component does not significantly affect the derived luminosity. For 9 l Ori C 
(COUP 809), the brightest X-ray source in the field, substantial discrepancies between best- 
fit model and data are left, which may be related to residual pile- up in the extracted area 
and/or uncertainties in the response for the annular extraction region. Since our modelling 
did not improve the result for COUP 809, we use the spectral parameters given by Getman 
et al. (2005a), but caution that the luminosity might be overestimated (see § A.l). 

The emission measures for the soft (EMi) and hard (EM2) components of the 2-T fits are 
compared in Figure 2. There is a general tendency of the WW stars to show EM 2 /EMi > 1, 
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and of the SW stars to show EM 2 /EM l < 1. Some SW stars - 6 2 Ori A, NUOri and to 
lesser extent Par 1772 - are very strongly dominated by the soft component. 

In addition to this plasma modelling procedure, we examined a simple model-independent 
spectral measure using the hardness ratio. The median values of HR for the three sub- 
samples are: (HR) M a = 0.37, (i?it!)ww = 0.29, and (HR) S w = 0.18. Nonparametric 
two-sample tests indicate that the difference in HR distributions for SW stars and MA stars 
is significant at the P < 1% confidence level, while WW and MA stars are statistically 
indist inguishable . 

4. Variability 

4.1. Statistical tests 

Variability is a characteristic property of magnetically active low-mass pre-main se- 
quence stars (Montmerle et al. 1983). Strong variations of their X-ray output are observed 
on short timescales (minutes to hours), while longterm variability is less pronounced and 
may be the result of catching the star during different phases of short-term events. Studies 
of X-ray variability on short timescales in hot wind-driven stars are practially absent in the 
literature. But comparing different sets of observations for a given star has provided little or 
no evidence for long-term variations (Berghofer et al. 1997), leading to the conviction that 
hot stars show no X-ray variability. In this vein, variability can be an important diagnostic 
for discriminating the origin of the X-ray emission in wind shocks vs. magnetic activity. On 
the other hand, to our knowledge no study has been performed where the variability of hot 
and cool stars was directly compared in a consistent manner. This is the purpose of this 
section. 

The observable variable phenomena include flares of different morphology, rotational 
modulation, and irregular variations on short and long timescales. For subsamples of COUP 
stars detailed variability studies are presented by Favata et al. (2005); Flaccomio et al. (2005); 
Wolk et al. (2005). Classification of a given event into any of the above mentioned categories 
of variability is not straightforward. Especially the definition of a flare, generally viewed as 
an event with rapid rise and spectral hardening followed by decay and spectral softening, 
escapes consensus on its quantitative description throughout the literature. 

We apply here a method described by Wolk et al. (2005) for analysis of flare-like vari- 
ability in COUP MA sources. It is based on Maximum Likelihood Blocks (MLB), a mod- 
ification of Bayesian Blocks developed by Scargle (1998), in which the X-ray lightcurve is 
divided into contiguous segments of constant signal. The segmentation algorithm has two 
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free paramenters: the minimum number of counts accepted in a given segment and the sig- 
nificance, defined as the probability that a constant lightcurve is spuriously segmented due 
to random fluctuations. In this study of OBA COUP stars, we set a minimum segment of 
20 counts and a significance for change points of 95%. 

The results of the variability tests are summarized in Table 4. Stars are sorted in order 
of decreasing effective temperature. Columns 4 and 5 flag stars variable according to a simple 
Kolmogorov-Smirnov (KS) test and the MLB test. The KS-test shows 12 out of 16 targets 
are variable at the 95 % confidence level. The MLB test indicates 14 of our targets can not 
be described with a single block and thus are likely variable. 

The ML blocks are also used to define a 'characteristic level' and flares, represented by 
periods of 'very elevated' count rate and rapid flux change as quantified by the derivative of 
the lightcurve (Wolk et al. 2005). A total of 15 flares are identified, distributed among the 
20 targets as summarized in column 6 of Table 4. Columns 7 and 8 give the fraction of the 
time spent in the 'characteristic level' (£ c har/£totai) and the ratio between the count rate of 
the highest segment and the 'characteristic level' (-R pe ak/-Rchar)- In Figure 3, the amplitude 
of the lightcurve variations, R max /R m i n is compared to that of COUP low-mass MA stars. 
Stars for which the maximum count rate corresponds to a flare are distinguished by filled 
plotting symbols. The amplitudes clearly indicate that SW stars are less strongly variable 
than MA stars. 

The sensitivity to the detection of variability depends on the signal-to-noise. Fig. 4 
shows the fraction of variable stars for different count levels. The curves represent COUP 
MA stars. The solid line is for the KS-test and the dashed line for the MLB test. The total 
number of stars in each bin is indicated on top of the graph. These data demonstrate that 
nearly all X-ray-bright low-mass stars are found to be variable. An equivalent study for the 
SW and WW stars is prohibited by the small number of objects in our sample. Instead, 
we place the SW and WW stars on the MLB and KS curves for MA stars, at the position 
corresponding to their observed counts. From this comparison, we find that 96% of WW 
and SW stars are expected to be variable. For the WW, this is consistent with the observed 
fraction of 100% variable stars. The fraction of variable WW stars is thus similar to that 
of MA stars. But for the SW stars, the expected variable fraction is somewhat higher than 
the observed fraction of 78% variable stars (see Table 4). A deficiency of SW variability is 
also suggested by the fact that only 4 flares appear in 9 SW stars compared to 11 flares in 
7 detected WW stars. 
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4.2. Flares 

The lightcurves of all detected targets are shown in Fig. 5 together with the time evo- 
lution of the hardness ratio evaluated in each segment of the MLB test. Signatures of 
heating during phases of increased intensity are obvious support for the interpretation of 
these episodes as magnetic reconnection flares. In the following, the characteristics of the 
flares on SW and WW stars are compared to the flares on solar-mass COUP stars presented 
by Wolk et al. (2005) which uses the same MLB and flare characterization methodology. 
Since only 4 flares were observed on SW stars, we did not compute statistics for this group 
alone, but rather examine how the derived values for WW stars change if we combine the 
sample with the SW stars. Spectral analysis is based on two-temperature plasma fits of flare 
segments (see § 3) for the 12 (9) flares in the SW+WW (WW) samples where >50 counts 
were collected. The results of this comparison are summarized in Table 5. 

A histogram showing the distribution of flare durations for the WW stars is given in 
Figure 6. Although the median of the flare duration on WW stars is only about 1/3 of 
the value for solar analogs, the distribution seems to be roughly in agreement with the 
0.1 x e (-*/ 100fcs ) relation that describes the distribution for solar-type stars. Following the 
arguments of Wolk et al. (2005), we derive a flare frequency of 1 flare every 710 ks for WW 
stars which is similar to the value of 1 flare every 640 ks obtained for ONC solar analogs. 
For SW stars, the flare frequency is three times lower at 1 flare every ~ 2200 ks. 

The median flare luminosity and energy of WW stars is similar to the respective values 
for solar-type stars, while in the total sample of SW plus WW stars both the flare luminosity 
and energy are increased by almost 1 dex. A possible interpretation is that the flares of 
SW stars are not drawn from the population of COUP MA stars, but are the result of a 
different emission mechanism. Since the sudden energy release observed in flares points at 
the presence of magnetic fields one possible explanation for the flares on SW stars is that 
they arise from magnetically mediated events in the strong winds. Note, however that the 
flares on SW stars are not particularly strong with respect to the characteristic level; their 
high luminosities are due to the overall higher emission level of SW stars with respect to the 
WW and solar-type stars. The different quiescent luminosity of the SW stars and WW stars 
may give rise to a bias, because smaller flares are more difficult to detect on brighter stars. 
Indeed, if we examine the individual events, the strongest flares on SW stars are comparable 
to the strongest flares on WW stars and solar-type stars, at odds with the above mentioned 
interpretation. A thorough quantitative comparison of COUP spectra/ variability and the 
MCWS model can be obtained with help of simulations, but it is beyond the scope of this 
paper and will be addressed in a future study. 
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4.3. Rotational Modulation 

Rotational modulation of the X-ray signal is expected if the X-ray emission is concen- 
trated in long-lived inhomogeneous magnetic field structures, such as active regions on a MA 
star or the magnetosphere in a magnetically confined wind shock (MCWS). However, since 
the X-ray emitting plasma is optically thin, rotational modulation may exist only as a result 
of asymmetries along the line of sight: stellar occultation (provided the emitting volume is 
not much larger than the star), or non-spherically symmetric absorbing material such as the 
'cooling disk' of the MCWS model. 

In contrast to the optical band where rotational patterns related to active regions are 
readily observed in many late-type stars, reports on rotational modulation in the X-ray band 
have been scarce due to inadequate time coverage of past observations. The rotation periods 
of active late-type stars are on the order of ~ 1 — 10 days, accessible to the 13.2 days spanned 
by the COUP exposure but longer than previous X-ray exposures of pre-main sequence stars. 

Flaccomio et al. (2005) detect rotational modulation in ~ 10 % of a subgroup of > 200 
COUP stars examined systematically, which definitely is the biggest sample studied in this 
respect to date. They compute the Lomb Normalized Periodogram (Scargle 1982) on binned 
COUP lightcurves. The peak frequencies in the periodogram are then identified and a False 
Alarm Probability (FAP) is computed using Monte Carlo techniques allowing for correlated 
noise. The analysis is made both on the original data and on lightcurves trimmed of bright 
flares. We now apply these methods to the ONC early-type star sample. 

We examine the COUP lightcurves of the three stars in the SW and WW sample which 
have photometrically measured rotation periods: 9 1 OriC, JW660, and JW165. Two addi- 
tional stars, JW831 and JW531, were surveyed for starspot induced photometric variability 
by Choi & Herbst (1996), but no period was found. At the present stage we discuss only 
stars with known optical period, because the search for X-ray periods is complicated by 
flares and other types of variability. A systematic analysis of the whole COUP sample with 
known optical period (Flaccomio et al. 2005) has shown that sometimes the X-ray period 
does not coincide with the optical one, suggesting that the complex X-ray lightcurves can 
result in spurious periods. Our results for the three stars analysed here are summarized 
in Table 6, and phase-folded lightcurves and hardness ratios can be found in Fig. 7. The 
median count rate and the amplitude (columns 5 and 6 of Table 6) were computed from the 
folded lightcurve. The uncertainty of the amplitude derives from the scatter of the data near 
maximum and minimum, and thus takes into account the intensity changes from one cycle 
to the other. 
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4.3.1. 9 1 On C = COUP 809 

An optical period of 15.43 d was reported by Stahl et al. (1993), and also evidenced in 
a set of 10 ROSAT X-ray observations separated by ~ 2 days each, such spanning a little 
more than one cycle (Gagne et al. 1997). The same variability pattern is recovered in the 
COUP lightcurve (Fig. 7). The X-ray lightcurve is almost exactly in phase with the Ha 
ephemeris measured 10 years earlier, and the X-ray period agrees remarkably well with the 
optical period. The X-ray hardness ratio varies out of phase with the intensity by roughly 
half a cycle. While there is a concern that incomplete correction for photon pile- up in the 
ACIS detector may be responsible for the variation of the hardness ratio, we believe this is 
not occurring. If this were the case, one would expect a correlation with the lightcurve and 
not an anti-correlation because the increased pile-up during higher intensities would shift 
photons into the high-energy portion of the spectrum. 

Astrophysical discussion of this source in the context of detailed evidence for an inclined 
dipole field in 9 l Ori C and magnetically trapped wind plasma (e.g. Babel & Montmerle 
1997a; Donati et al. 2002; Schulz et al. 2003; Smith & Fullerton 2005) will be the subject of 
a forthcoming study. 

4.3.2. JW 660 = COUP 1116 

A rotation period of 6.15 d was reported by Choi & Herbst (1996) and Herbst et al. 
(2000) for this B2.5 star, although the cause of a spot pattern on an early-B star is not clear. 
We find two episodes of enhanced count rates in the COUP lightcurve, which may define an 
X-ray period of 6.67 d (Fig. 7). The amplitude of the folded lightcurve is poorly defined due 
to the different peak intensity in the two consecutive cycles. It is also possible that the two 
features are long-duration flares unrelated to any rotating structures. The hardness ratio 
varies in phase with the brightness of the source. Spectral hardening during phases of higher 
intensity is a typical signature of magnetic activity and is opposite to the pattern seen in 
^OriC. 

More detailed information on the run of temperature and emission measure is obtained 
from time-resolved spectroscopy. The spectrum of JW660 can be described adequately 
by a two-temperature thermal model with comparatively high temperature at all phases 
of the presumed rotation cycle. These fits require an absorbing column somewhat higher 
than the value inferred from the optical extinction. During the maximum of the cycle the 
higher temperature is slightly enhanced and the emission measure shifted towards the hot 
component, in agreement with the evolution of the hardness ratio. A feature from highly 
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ionized iron at ~ 6.7 keV is pronounced during the phase of maximum brightness. Its 
appearance results from improved statistics rather than a true change in the spectral shape. 

4.3.3. JW 165 = COUP 113 

A photometric period of 5.77 d was reported by Herbst et al. (2000) for this A7 star, 
albeit based on a single observing season. The period search on the COUP data results in 
a period of 4.6 d, but the folded X-ray lightcurve is rather noisy (Fig. 7). Spectral fitting 
suggests a soft X-ray absorption inconsistent with the visual extinction. 

5. Multiple Systems 

About half of our sample has been carefully examined for multiplicity using traditional 
spectroscopic and newer speckle and adaptive optics imaging techniques, as indicated in 
columns 7 — 9 of Table 1. High-multiplicity visual systems are common among the ONC OB 
stars, and we summarize known components of our sample stars in Table 7. In all cases the 
primary falls into the SW category, and most of the secondaries seem to be WW stars rather 
than MA stars. However, the mass estimates are often based on photometry and have large 
uncertainties. We also caution that additional, faint lower-mass companions may be present 
which could contribute to the X-ray emission. Recall that 1" separation corresponds to a 
projected distance of 440 AU. Five well-characterized systems are presented in Appendix A 
where we discuss the characteristics of their X-ray emission with the aim to identify the 
X-ray emitting component. 

Only one of these multiple systems is resolved in the COUP: the quintet 9 1 Ori B consists 
of two X-ray sources (Fig. 8, see ApppendixB for details). ^OriBEast (COUP 778) emits 
logL t ,c = 30.1 erg s _1 with an unusually soft spectrum and no strong variability. This source 
can be plausibly attributed to wind shocks of the B2.5 primary. In contrast, 9 1 Ori B West 
(COUP 766) has all the characteristics of a MA lower-mass member of the Orion population 
with several flares. It is not clear which of the components of the 9 l Ori B West binary is the 
principal X-ray emitter. 

6. Stars with evidence for disks 

Our sample includes four stars listed in the catalog of Herbig Ae/Be stars by The et al. 
(1994) in their Table 5 labelled 'non emission line early-type shell stars and young stellar 
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candidates'. Stars in this table lack emission line signatures of accretion, thus do not satisfy 
the classical definition of HAeBe stars as emission line objects, but IR photometric excesses 
were seen with IRAS. A dust shell is identified around one additional star of our sample 
(Tovmassian et al. 1997, see flag for IRAS excess stars in the last column of Table 1). The 
IRAS measurements indicate that a substantial amount of the luminosity of these stars 
is irradiated in the far-IR, probably related to dust emission at 100 — 1000 AU from the 
star. Par 1772 and V361 Ori are discussed by Manoj et al. (2002) who estimate dust masses 
between those of typical HAeBe stars and of Vega-like stars. The absence of emission features 
in the spectrum of Par 1772 and JW153 was confirmed by Bohm & Catala (1995), thereby 
showing that the accretion phase has ended in these objects. 

Figure 11 shows the near-IR color-color diagrams for our sample. For those stars de- 
tected in the COUP the JHKL magnitudes are reported by Getman et al. (2005a) while for 
the nondetections we compiled IR photometry from the 2 MASS database and the literature 
(Hillenbrand et al. 1998; Carpenter et al. 2001; Muench et al. 2002). No strong indications 
for excess emission above the photosphere is seen in the Jffi^-diagram, but several stars 
have L-band excesses. A comparison to the intrinsic IR colors of a late-type star shows that 
the redder color expected from the presence of a hidden cool star can not account for the 
observed excesses. We conclude that the near-IR excesses must be ascribed to circumstellar 
material. We note that the five IRAS excess stars do not show stronger K — L excess than 
the remaining sample. This suggests they have evolved disks where the inner, hotter part 
has been cleared. 

We find no link between X-ray properties and the presence or absence of an IR excess. 
For example, a wide spread —8 < log (L x /Lboi) < —5 is seen among the stars with evidence 
for disks, which does not favor a common origin of their X-rays (Figure 12). One of the five 
IRAS excess stars, JW 531, remains undetected in the COUP. 

7. Discussion 

7.1. COUP results in perspective 

The COUP is the most sensitive X-ray image ever obtained of a rich young (~ 1 Myr 
old) stellar cluster with an initial mass function extending up to ~ 45 M , providing the 
opportunity to study the X-ray emission of a coeval sample of early-type stars on a long 
continuous time baseline. The X-ray emission mechanisms of early-type stars are highly 
disputed (§1). In the hotter stars, our SW sample, the stellar winds most likely play a 
crucial role. About half of the SW sample is dominated by soft X-rays that are likely 
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generated by a myriad of small shocks produced by instabilities in their line-driven winds 
(Lucy & White 1980; Owocki k Cohen 1999). The other half, mostly the Bl to B3 type 
stars, the cooler stars in the SW sample, also show substantial hard emission, unexplained 
within the wind-shock model. This emission may indicate that magnetic fields are involved. 

Direct observational evidence for the existence of magnetic fields on early-type stars is 
scarce and mostly restricted to the class of CP stars (Landstreet 1992). On the other hand, 
magnetic fields are often held responsible for imposing the cyclic variability of wind proper- 
ties observed in the UV (e.g. Kaper et al. 1999; de Jong et al. 2001) and for the non-thermal 
radio emission seen in some O stars (Abbott et al. 1984; Van Loo et al. 2004). While tradi- 
tionally magnetic fields on hot stars are believed to be fossil remnants of the star formation 
process, numerical modelling has shown that dynamos may operate in their convective cores 
(Charbonneau & MacGregor 2001). However, the theory has difficulties transporting mag- 
netic flux throughout the thick radiative layer to the stellar surface (MacGregor & Cassinelli 

2003) . Hybrid models including winds disturbed in large-scale magnetospheres have been 
developed to explain both radio and X-ray emission from early-type stars (Trigilio et al. 

2004) . 

These models have two implications for the mechanism of X-ray production. First, X- 
ray emission may be produced as a result of magnetic reconnection in the current sheet that 
forms in the equatorial plane where gas pressure opens the field lines. This model explains 
the observed non-thermal radio emission of magnetic CP stars (Linsky et al. 1992), and 
Usov & Melrose (1992) have estimated that the expected X-ray luminosities are logL x ~ 
31 — 34 erg s -1 . Therefore, gas heating by magnetic reconnection or other instabilities at 
the equatorial current sheet from a dipolar field of an early-type star can account for X-ray 
emission from radiative stars; i.e., both SW and WW stars. Second, in the magnetically 
confined wind shock (MCWS) model (Babel & Montmerle 1997a), the major effect of the 
magnetosphere is to confine and channel the wind to the equatorial plane where a quasi- 
stationary large-scale shock forms. The cooling of shocked material results in the formation 
of a 'cooling disk', partially absorbing the emitted X-rays and periodically modulating the 
observed X-rays. This model can explain the high temperatures and unusual X-ray line 
shapes seen in 9 1 Ori C (Schulz et al. 2003; Gagne et al. 2005). 

In less luminous WW stars, the wind velocities are smaller than in SW stars. Production 
of X-ray emission in small-scale wind shocks becomes questionable although the MCWS 
model may still apply; indeed, this model was first devised for the magnetic Ap star, IQ Aur 
(Babel & Montmerle 1997a). Perhaps the most common explanation for the X-ray emission 
from WW stars - applicable to the weaker SW stars as well - invokes unknown or unresolved 
late- type companions (separations smaller than ~100AU). Such objects would by their 
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nature be magnetically active (MA), and produce relatively hard, highly variable X-ray 
emission with typical luminosity between logL x ~ 28 — 31 ergs -1 . 

Distinguishing between the different possible emission scenarios based on observational 
facts is not straightforward. In this section, we use the results from the COUP to test our a 
priori division of the early-type stars in the ONC into SW and WW stars, and we provide 
a tentative classification of all targets. This classification scheme involves three criteria - 
X-ray luminosity, spectrum, and variability - and is summarized in Table 8. 

7.2. Classification of ONC early-type stars X-ray emission 

We consider first Fig. 12 which plots the L x /L ho \ ratio against the effective tempera- 
ture for the SW and WW stars in the COUP. For both groups of stars, a wide spread of 
L x /Lboi ratios is seen without dependence on T e g. The WW stars show a very wide spread 
around log (L x /Lb i) ~ —5 to —6. The SW stars do not cluster near the empirical relation 
log (L x /Lboi) — 1CT 7 reported in earlier observations of X-ray emission from hot stars and 
interpreted in terms of the standard wind-shock model (Sciortino et al. 1990; Berghofer et al. 
1997). They scatter by more than two dex around this level. 

Nonetheless, we use this L x — Lboi scaling relation to help classify and understand the 
X-ray emission of ONC early-type stars. We compare the observed luminosities L t;C to the 
'predicted' luminosity L XjP = 10~ 7 Lboi for stars of the SW class. For the WW class, we 
use a 'predicted' luminosity L XjP = 1CT 4 L bo i which is characteristic of lower-mass MA stars 
(Preibisch et al. 2005). These levels are shown as dotted lines in Figure 12. Since MA 
companions are likely to have L bo i lower than the primary, this L x p is actually a lower limit. 
In column 4 of Table 8, we distinguish the following situations: 

If L tjC /L XjP < 10, the emission mechanism defined above for the SW and WW stars respec- 
tively is considered to be the dominating mechanism. We denote this mechanism in 
column 4 by "W" for wind-shock and "T" for T Tauri-like magnetic activity. Note 
that, in the latter case, as an alternative to a MA late-type companion, the emission 
might be due to a magnetic reconnection process in a wind. 

If L tjC /L XiP > 10, a different or additional mechanism is considered to be present. We 
denote this situation by a ">" sign in column 4. 

In column 5 of Table 8, we summarize the result of the X-ray spectral analysis (§3) 
for all detected stars. We consider the X-ray spectrum characteristic for small-scale wind 
shocks if the soft component (with T ~ 2.5 MK) of the two-temperature fits dominates the 
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emission measure (see Figure 2), or if a hard component is altogether absent. Such sources 
are denoted by "W" . We see a clear pattern that such soft-dominated spectra are seen only 
in some SW stars and not in any WW stars. 

Column 6 tabulates results of the variability analysis (§4) based on the COUP light 
curves. Constant or weakly variable light curves characteristic of small-scale wind shocks are 
denoted "W". Flare-like high-amplitude variability indicates magnetic reconnection ("T"), 
which might be either from lower-mass MA T Tauri companions or from large-scale shock 
events in a magnetically confined wind. Smooth lightcurves ("S") indicate non-flare-like 
variability such as periodic modulations with moderate amplitudes (factors ~ 2 — 3). 

In column 7, we categorize the observational indications for the presence of circumstel- 
lar material (§6). Recall that no emission-line stars (= HAeBe stars) have been identified 
in the ONC, but some evidence for disks is seen in the near-IR color-color diagrams or in 
far-IR fluxes. Finally, column 8 collects observational evidence for magnetic fields includ- 
ing rotational modulations (RM), non-thermal cm radio emission (NT), chemically peculiar 
characteristics (CP), or direct indicators such as by Zeeman or polarization measurements 
(B). 

Some immediate conclusions may be drawn from Table 8 and the underlying findings 
presented throughout this paper: 

1. There is a clear distinction between the higher- mass SW and all WW stars in terms 
of their X-ray properties. This justifies a posteriori our adopted division between SW 
and WW classification. However, the location of the transition between the two classes 
remains ill-defined and may be earlier than the B4 boundary adopted in § 2.1. Some 
BO— B3 stars have luminosities consistent with X-ray production in the wind, but have 
spectra and variability characteristics of magnetic activity. 

2. In the more massive SW stars with spectral types BO— 07, small-scale wind shocks 
predicted by the original model of Lucy & White (1980) appear to be the dominant X- 
ray emission mechanism. But even here, magnetic processes seem to be present; e.g., to 
explain the powerful L x ~ 10 32 erg s" 1 COUP flare in 9 1 Ori A (Figure 5), and a similar 
flare in # 2 OriA seen in the earlier Chandra study of Feigelson et al. (2002). If one 
includes the rotational modulation of 9 1 Ori C and the X-ray independent indicators 
of magnetic fields, all three of the most massive Trapezium stars - 6' 1 OriC, # 2 OriA 
and 8 l Ori A - show hybrid wind/magnetic characteristics; see Schulz et al. (2003) for 
a discussion of these hybrid systems in the Orion Nebula hot stars. These stars are 
promising candidates for the MCWS model. 

3. The less massive SW stars similarly show hybrid characteristics combining winds and 
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magnetic activity. Here the X-ray luminosities are lower, L x ~ 29 — 31 erg s -1 , as 
expected for lower-mass T Tauri companions, and evidence for magnetic fields is found 
in some but not all objects. Interestingly, most SW stars with emission above the 
canonical line of log (L x /L bo i)= —7 are known to have lower-mass companions that are 
unresolved in X-rays. 

4. In the WW class, the X-ray emission is either undetected with very sensitive limits 
(L x /Lboi < 10~ 7 ), or is detected with properties associated with magnetic activity from 
a late-type companion. 

5. Objects with circumstellar disk indicators show no evidence for a common X-ray emis- 
sion mechanism. We conclude that disks do not play a significant role in the X-ray 
emission of the early-type ONC stars. 

We will now discuss the sample stars in more detail, starting with the cooler intermediate- 
mass WW stars and proceeding to the SW stars. 

7.3. WW stars 

The situation for the ONC WW stars (i.e., spectral types B5 and later) seems fairly 
straightforward (Table 8). Out of 11 stars, 4 are not detected with upper limits several 
orders of magnitude below the activity level and clearly lower than the L x of all detected 
WW stars (see Fig. 12). A bimodal distribution of X-ray luminosities for B- and A-type 
stars has been suggested in other studies (Preibisch & Zinnecker 2001; Feigelson et al. 2002; 
Stelzer et al. 2003), explained by zero emission by the intermediate-mass stars and emission 
by lower- mass MA companions when present. Our data can not distinguish between two 
separate populations and a single population with a very wide dispersion. 

The X-ray luminosities of the majority of the detected WW stars range between log L x [erg s _1 ] ~ 
29 — 31 erg s" 1 , typical of ONC late-type MA stars (Feigelson et al. 2005; Preibisch et al. 
2005), and therefore consistent with the unresolved binary companion hypothesis. Our study 
of their variability in general and flaring in particular, and their spectral hardness support the 
identification of WW X-ray emission with late-type MA companions. Not enough is known 
about 9 1 Ori BW to assign an X-ray emission mechanism, although the flare-like variability 
suggests the possible presence of an MA companion. Remarkably, none of the WW stars in 
the COUP sample is a known multiple; the multiple systems discussed in Appendix A are 
all SW stars. However, this is most likely due to a lack of systematic binarity searches in 
the WW population. 
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The WW star 6 1 OriE (B5) has unique properties in the COUP data. It has a much 
higher time-averaged X-ray luminosity than the other WW stars, logL t c = 32.4 erg s -1 and 
log(L t c /L bo i) = —3.6 (Table 3). Its log (L x /L bol ) level is several orders of magnitude above 
any standard wind-shock prediction. This star also shows evidence for non-thermal radio 
emission (Felli et al. 1993), and Schulz et al. (2003) found plasma temperatures ranging from 
4 to 47 MK in a study of Chandra grating observations. Therefore, it may have an extended 
magnetosphere, able to confine its weak wind within the MCWS framework. ^ 1 OriE has 
also shown a flare, which - if attributed to a late-type MA companion - would be one of 
the most powerful flares seen in any T Tauri star (Grosso et al. 2004; Favata et al. 2005). 
Altogether, the COUP properties of 6 1 Ori E are not easily explained in a coherent fashion. 

The X-ray properties of the A7 star JW 165 are also worth noting. Its smooth lightcurve 
is consistent with both rotational modulation or slow flares. Its late spectral type implies 
that its wind must be very weak. For comparison, the calculations by Babel & Montmerle 
(1997a) for IQ Aur, a hotter AO star, give M ~ 1O~ 1O M yr _1 . JW165 may have an unseen 
magnetosphere, although in this case its mass-loss rate should be M > 1O~ 12 M yr" 1 since 
it is not known to be CP. Perhaps more likely, it has a MA T Tauri companion responsible 
for the emission and peculiar variability observed during the COUP exposure. 

7.4. SW stars 

Based on the emission diagnostics in column 4 of Table 8, the SW stars can be subdivided 
in two groups: (i) those for which the X-ray luminosity is consistent with small-scale wind 
shock emission, and (ii) those which are overluminous by one or two orders of magnitude 
with respect to the corresponding predictions. With the exception of JW 660, all stars in 
the B0.5-B3 range fall in the "W" category according to their X-ray luminosity. However, 
a MA companion is suggested for # 2 OriB. Here the log (L x /L bo i) level is compatible with 
the prediction for wind-shock emission, but in absolute terms L x is low enough to originate 
from a late-type T Tauri star. 

JW 660 (B3) is one of three stars showing evidence for rotational modulation in the 
COUP. The time-evolution of its hardness ratio shows a correlation with the X-ray brightness, 
opposite to the anticorrelation found in 9 l Ori C where the MCWS mechanism appears to 
dominate. Rotational modulation of the X-ray emission can result from the eclipse of a 
large magnetic loop by the star (Stelzer et al. 1999), and the joint modulation of intensity 
and temperature is typical for magnetic activity phenomena. However, we can not exclude 
applying the MCWS scenario to JW660, as this model does not make predictions on the 
variation of hardness with phase. It would be very valuable to search JW 660 for a large-scale 
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magnetic field. 

The remaining stars (BO-07) are the most massive in the ONC. They are all very lu- 
minous with logL t c > 32 erg s -1 , about one order of magnitude or more above the value 
"predicted" for small-scale wind shocks. Next we discuss these stars in turn. 

The most extreme case is obviously ^ 1 OriC with logL tc ~ 33.3 erg s _1 . ROSAT 
discovered an X-ray modulation with the same period as the optical period (P = 15. 4d; 
Gagne et al. 1997). To explain its extreme X-ray luminosity and rotational modulation, Babel 
& Montmerle (1997b) applied a MCWS model originally devised for the X-ray emission of 
the Ap star IQ Aur (Babel & Montmerle 1997a). For 6 l OriC, the MCWS model postulates 
a magnetic field of order 400 G buried in the wind and confining its dense, low-velocity 
layers. Because the outer layers of the wind were predicted to be open and sitting above the 
magnetically confined layers, such a magnetic field was difficult to detect. However, this was 
achieved by Donati et al. (2002) who found the predicted field via a rotationally modulated 
Zeeman effect. The COUP data enables for the first time a continuous study of X-ray spectral 
changes along the rotation cycle of 9 1 Ori C. The anti-correlation between source luminosity 
and hardness is unusual and detailed modelling is needed. At present, we speculate that 
the explanation for this anti-correlation results from particular viewing angles of the X-ray 
absorbing 'cooling disk' predicted by the MCWS model. In a recent study based on high- 
resolution spectroscopic X-ray observations of 9 1 Ori C Gagne et al. (2005) reported a small 
increase of the column density during phases when the disk is viewed edge-on, corresponding 
to the X-ray minimum. This is consistent with the increase in spectral hardness we observed 
during the COUP. 

The second-most hot star, 6 2 Ori A (09.5), can be accommodated in the class of small- 
scale wind shock emitters according to its characteristics shown during the COUP, but 
dramatic flares were seen in the earlier Chandra observation described by Feigelson et al. 
(2002). These flares can be caused by magnetic reconnection in the wind (which implies that 
this star be magnetic), or by an unusually active unseen MA companion. 

9 l Ori A similarly shows hybrid characteristics. This star shows evidence for magnetic 
fields from its non-thermal radio emission (Felli et al. 1993), and an X-ray spectroscopic 
study shows plasma with temperatures up to 43 MK (Schulz et al. 2003). With its flare-like 
behavior, it reaches X-ray luminosities unusually high for a low-mass MA companion. We 
speculate that its magnetic field is too weak for the wind to be confined, but that strong 
reconnection events take place along the wind. If the star rotates, the wind flow will, like 
the solar wind, take the shape of an Archmedean spiral past the Alfven radius. In such a 
configuration, "corotating interaction regions" (CIRs) analogous to those in the solar wind 
(e.g., Kissmann et al. 2004), or in the wind of the young pole-on B star AB Aur (Catala 
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et al. 1999), may exist. These CIRs could trigger reconnection events, such as proposed by 
Usov & Melrose (1992), to which the flare-like events observed on 9 l Ori A by COUP might 
be related. A direct search for magnetic fields on the stellar surface would be valuable. We 
recall, however, that 8 1 Ori A has two known lower- mass companions, and possibly others 
still undetected. The observed X-rays could arise from one extraordinarily MA companion, 
similar perhaps to LkHa 312 seen elsewhere in the Orion molecular cloud (Grosso et al. 
2004). 

8. Conclusions 

The COUP observation has permitted an in-depth study of the X-ray emission of the 
early-type stars of the ONC with the unprecedented exposure covering most of 13.2 con- 
tiguous days. We examine a sample of 20 coeval O, B and A stars. We address the X-ray 
properties and underlying astrophysical processes of two subclasses: 9 hotter members with 
spectral types 07— B3, designated 'strong wind' or SW stars, and 11 cooler members with 
spectral types B5— A9, designated 'weak wind' or WW stars. We also investigate the effects 
of multiplicity, disks, and independent evidence for stellar magnetic fields. 

Addressing the mystery of X-ray emission from intermediate-mass late-B and A stars, we 
find an extraordinary range in luminosities ranging from nondetections with logL x < 27.6 
erg s _1 comparable to the quiet Sun to violently flaring sources up to logLx ^ 31 erg 
s _1 . Based on flaring lightcurves and hard spectra, the detected WW stars are very likely 
associated with unseen late-type companions. The case of 9 1 Ori E (B5) is unusual with 
strong constant X-ray luminosity around Lx — 2 x 10 32 erg s _1 and flares of comparable 
intensity. Exhibiting non-thermal radio emission, this may be a case for MCWS emission in 
a mid-B star. 

The situation with the Trapezium SW OB stars is rich and intriguing. The X-ray 
emission gives little support for the long-standing claim that L x /Lb i is constant around 
~ 10~ 7 for OB stars. While a rough trend between L x and Lboi is seen, three orders of 
magnitude scatter in L x can be present at a given Lboi 3 - We note that this result is limited 
by the Orion population to stars of spectral type 07 and later, and we can not say whether 
a constant L x / L bo i is present for the most luminous early-0 stars. 

The X-ray properties from the majority of the SW stars exhibit hybrid characteristics. 



3 We suggest that the previous evidence, particularly the R OS A T All- Sky Survey of OB stars (Berghofer 
et al. 1997), suffered a distance-related selection bias and missed the lower luminosity stars. 
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Only 9 1 OriD (B0.5) and NUOri (Bl), and perhaps the faint Par 1772 (B2) and 6 1 OriBEast 
(B3), exclusively show the roughly constant, soft spectrum emission attributable to the tra- 
ditional Lucy- White model of small-scale wind shocks. Most show combinations of soft wind 
emission, hard-spectrum short-lived flares, and/or rotational modulation. If the flares of 
9 2 Ori A from a previous Chandra observation are included, all three of the hottest Trapez- 
ium stars - ^OriC (07), 6> 2 0riA (09.5) and ^OriA (BO) - show these hybrid char- 
acteristics. Together with independent indicators of magnetic fields, the COUP evidence 
strongly supports the MCWS model where the harder and variable X-rays are produced by 
magnetically-mediated large-scale shocks. 

The eleven cooler WW stars lack both strong ultraviolet emission to drive powerful winds 
and outer convection zones to drive a solar-like dynamo. They are thus predicted to be X-ray- 
quiet. Yet seven of these stars are detected in the COUP image and together they exhibit a 
wide range of X-ray luminosities from logL x < 27.6 erg s _1 comparable to the contemporary 
Sun to flares up to logL x ~ 31 erg s^ 1 . Based on their flaring lightcurves and hard spectra, 
very similar to COUP properties of lower-mass stars in the ONC, we believe that the WW 
X-ray emission is most likely produced by late-type companions. In one case, 6 1 OriB, the 
high-multiplicity system is resolved into two X-ray sources; one component seems to be a soft- 
spectrum wind source and the other one a flaring MA lower-mass companion. The COUP 
study uncovered a similar multiple system, heavily obscured by the Orion Molecular Cloud: 
the hot Becklin-Neugebauer Object has a faint X-ray source that is nearly overwhelmed by 
a much brighter, previously unknown low- mass companion (Grosso et al. 2005). 

Our boundary at spectral type B4 between SW and WW stars is probably not a strict 
divide between wind and companion X-ray emission in all stars. The case of 9 1 Ori E at B5 
is unusual, showing strong constant X-ray emission around L x ~ 2 10 32 ergs _1 and flares of 
comparable intensity. It may represent a case for MCWS emission in a mid-B star. 

We emerge from the present study with an increasing recognition of the important role 
played by magnetic fields in massive stellar systems. In many cases, the X-ray properties 
no longer follow the simple ideas of small-scale shocks from unstable radiatively accelerated 
winds. Violent reconnection events, either in the winds of the massive primary or from 
MA orbiting late-type companions, are frequently present. Rotational modulation of the 
X-ray signal is not infrequent. The high fraction of OB stars which are suggested from 
X-ray emission to have magnetically confined winds is much higher than the few percent 
of B stars possessing optically detected magnetic surfaces (Landstreet 2003). X-rays thus 
provide a powerful diagnostic for OB magnetospheres and their interactions with radiatively 
accelerated winds. The evidence provided here concerning the prevalence of strong magnetic 
fields in the hotter 07— BO stars, where optical detection techniques are difficult, may be 
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particularly important. 

It is now necessary to search for unresolved late-type companions down to distances ~ 
10 AU, as is currently possible with high-angular resolution optical and IR methods. On the 
other hand, our study must encourage renewed efforts to search for magnetic fields in OB 
stars, in particular in the Orion Trapezium. 

A. Multiple stellar systems 

Here we present observational information on individual multiple systems in order of 
decreasing spectral type of the primary component. 

A.l. ^OriC - COUP 809 

6 1 OriC is a speckle binary with 0.037" separation (or 16 AU projected linear separa- 
tion), composed of the 45 M Q primary and a > 6 M secondary (Schertl et al. 2003). It is 
the most massive star in the Trapezium cluster, and also extraordinary in terms of its X-ray 
properties: highest X-ray luminosity among the COUP sources (logL tiC [erg/s] ~ 33.3), rota- 
tional modulation of the X-ray emission in accordance with the optical variations, hardness 
anti-correlated with the lightcurve along the 15.4 d cycle. These pecularities make it likely 
that the X-ray emission originates on the O star. 

Although we took carefully care of correcting for pile-up effects, the extracted spectrum 
of this extremely X-ray luminous source may still be affected by pile-up, so that the spec- 
tral parameters and the luminosity are possibly misestimated. Schulz et al. (2001) derived 
log L x [ergs -1 ] ~ 32.3 from the zeroth order signal of an ACIS/HETG observation. But 
Gagne et al. (2005) found logL x [ergs -1 ] ~ 33 from the same data set, in rough agreement 
with the COUP result. Furthermore, Flaccomio et al. (2003) found logL x [ergs -1 ] ~ 33.2, 
consistent with our COUP result, from a Chandra HRC observation, which is free of pile-up 
effects. 

A. 2. # 2 OriA - COUP 1232 

This system is a hierarchical triple with a 25 M 09.5 primary (Preibisch et al. 1999), 
a ~ 9 M A2 spectroscopic companion (Aikman et al. 1974; Abt et al. 1991; Morrell & 
Levato 1991) and a ~ 7 M visual companion with separation 0.38 "(Mason et al. 1998; 
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Preibisch et al. 1999). The X-ray emission is unusually strong \ogL tiC = 32.5 erg s _1 after 
correction for absorption. The COUP source suffers heavy pileup and only 1% of the events 
are extracted for analysis. 

Such high levels of activity are rarely seen in lower mass T Tauri stars and, in such cases, 
are always associated with high amplitude (L flarc /L cha r >> 10), hard spectrum (kT ~ 5 — 8 
keV) flares (e.g., Tsuboi et al. 1998; Grosso et al. 2004; Feigelson et al. 2005; Favata et al. 
2005). However, 6 2 Ori A is only weakly variable over the 13-day COUP observation: no flare 
is seen and the luminosity ratio between the highest segment and the characteristic level is 
-kpeak/^char = 1-8. In the earlier Chandra ACIS study of the ONC of Feigelson et al. (2002), 
8 2 Ori A exhibited multiple flares with peak logL t ~ 31.8 erg s _1 superposed on a relatively 
constant component. This was the most dramatic variation then known from an OB stellar 
system and can be attributed either to rapidly changing large-scale shocks in a magnetically 
confined stellar wind of the massive primary, or to unusually powerful magnetic reconnection 
flares in a lower mass companion. 

The spectrum is well- fit by a two-temperature plasma dominated by the soft component. 
The absorption column derived from the X-ray spectrum is compatible with the optical 
extinction. We thus believe the emission arises from the massive primary wind and not a 
MA secondary. Note that our individually constructed X-ray spectral fit differs from that 
obtained by Getman et al. (2005a) from a semi-automated procedure. They reported a much 
lower absorbing column density, and hence lower value of the absorption-corrected logL t c . 

A.3. 9 1 Ori A - COUP 745 

9 1 Ori A is at least a triple system with a 16 M B0. 5 primary (Schertl et al. 2003). The 
nearer companion has an eclipsing 65 day orbit; models of the optical lightcurve suggests a 
pre-main sequence ~ 3 M ~A0 star (Lloyd & Strickland 1999). The distant companion is a 
4 M late-B star with 0.2" separation (Petr et al. 1998; Weigelt et al. 1999; Preibisch et al. 
1999). The system is a highly variable radio source with spectral and polarization properties 
characteristic of non-thermal flaring (Felli et al. 1993). 

9 l Ori A is a strong X-ray source which is affected by photon pile-up in the ACIS detector. 
During the 13-day COUP observation, it exhibited a nearly-constant 'characteristic' level of 
log L tjC ~ 32.4 erg s -1 superposed with several flares with peak luminosities around 31.9—32.6 
erg s _1 . A range of plasmas with energies from 0.7 to >2 keV are present. 
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A.4. NU Ori - COUP 1468 

NU Ori is at least a hierarchical triple binary with a 14 M Bl primary, a ~ 3M Q 
spectroscopic companion with period of 8 days (Morrell & Levato 1991), and a ~ 1M 
companion with separation 0.47" (Preibisch et al. 1999). The COUP X-ray luminosity is 
logL t c = 31.2 erg s" 1 . This source was located 8' off-axis where the Chandra point spread 
function is degraded, so the 0.4" offset between the COUP and primary component positions 
can be reconciled with emission from any of the components. However, the L x /Lboi level 
is that of a typical wind X-ray source, and NU Ori showed no sign of variability during the 
13-day observation. Together, these characteristics suggest that the observed emission comes 
from the primary and can be explained within the standard scenario of stellar winds. 

A. 5. 9 1 Ori B - COUP 766+778 

This system is at least a quintet. Near-IR images of the system are provided by Simon 
et al. (1999) and by Schertl et al. (2003). The brightest member is the B2.5 component 9 l 
Ori B East with mass 6.3 M (Prato et al. 2002). It is an eclipsing binary (hence the variable 
star designation BM Ori). The secondary is a rapidly rotating pre-main sequence star with 
large IR excess, estimated spectral type of G2 (Popper & Pavec 1976; Antokhina et al. 1989; 
Vitrichenko & Plachinda 2000), and mass of ~ 2.7 M (Palla & Stahler 2001). 9 1 OriB West 
lies 0.94" to the southwest and is a resolved binary with 0.14" separation. The estimates for 
the masses of these two components range between 4 — 1.6 M and 3 — 0.7 M , respectively 
(Schertl et al. 2003; Preibisch et al. 1999). The faintest known member of the quintet has 
i^-band magnitude 10.5 (Simon et al. 1999) with estimated mass 0.2 M (Preibisch et al. 
1999), and lies ~ 0.6" northwest of 9 l OriB East. 

9 1 Ori B is the only visual multiple that is resolved into multiple COUP sources, although 
several dozen COUP doubles with separations < 3" are seen in the MA population (Getman 
et al. 2005b). The COUP image of the region, shown in Figure 8, shows two X-ray sources: 
9 1 OriB East (COUP 778) and West (COUP 766). About 1/3 of the events from each source 
is retrieved with the ACIS Extract procedure to reduce cross-contamination of the point 
spread functions. However, some residual contamination of the luminous flare in the brighter 
source 766 is seen in the lightcurve of the weaker source 778 (Figure 5). 

^ 1 OriBEast (COUP 778) with logL tiC = 30.1 erg s _1 has an unusually soft spectrum 
without strong variability. The emission is thus plausibly attributed to wind shocks of the 
B2.5 primary, rather than to reconnection flares in its G2-type companion. 

In contrast, 9 1 Ori B West (COUP 766) has all the characteristics of a magnetically active 
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lower-mass member of the Orion population (Wolk et al. 2005; Favata et al. 2005). It shows 
several flares including a very powerful one exhibiting rapid rise with a hard spectrum, peak 
luminosity around logL t ~ 31.5 erg s -1 , and slower decay over several hours with a softening 
spectrum. It is not clear which of the components of the 8 1 Ori B West binary is the principal 
X-ray emitter. 

B. New X-ray Detection of Par 1772 ( = COUP 349) 

The high spatial resolution of Chandra combined with the high sensitivity achieved 
during the extraordinary long exposure of the COUP has brought forth a new X-ray detected 
B1.5 star, Par 1772. It is sometimes considered a Herbig Be star due to mid-IR excess in 
the IRAS satellite data and optical reddening. However, it does not exhibit emission lines 
indicating accretion nor near-IR photometric excess emission indicating a warm inner disk, 
as seen in other Herbig Be stars (Manoj et al. 2002). It may have an annular disk. No 
companions were seen in speckle imaging (Preibisch et al. 1999). 

In X-rays, the star was reported at a level of logL s ~ 30.2 erg s" 1 from ROSAT HRI 
observations in 1991-92 (Caillault et al. 1994). However, the offset was reported to be 7" so it 
seems likely that the association was confused by the three COUP sources within ~ 10 "(see 
Fig. 10). The star was not detected in 1999-2000 Chandra ACIS or HRC observations 
(Feigelson et al. 2002; Flaccomio et al. 2003). During the COUP 2003 exposure, it appeared 
as a weak and roughly constant source with logL tiC = 28.9 erg -1 . Together with 8 2 Ori B, 
it has one of the lowest log (L x /L bol ) values (= —7.9) ever seen in a stellar X-ray source. 
We can not determine whether the emission arises from the massive primary or an unseen 
companion. 
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suggestions of the referee M. Gagne. COUP is supported by Chandra Guest Observer grand 
SAO GO3-4009A (E. Feigelson, PI). B.S., E.F., CM. and S.S. acknowledge financial support 
from the Italian MIUR FRIN program and an INAF program for the years 2002-04. E.D.F. 
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Table I. Stellar parameters of early-type stars in the COUP field of the ONC 



ID 


V 


A v 


SpTy logT cfr 


log L bo i 




Multiplicity 






mag 


mag 




K 


L Q 


SB? a 


VB? b Optical? c 


Disk? d 




'Strong- Wind' Sample 


6 1 Ori C 


5.140 


1.74 


07 


4.603 


5.40 




X 




2 Ori A 


5.070 


1.12 


09.5 


4.543 


5.03 


X 


X 




6 l Ori A 


6.730 


1.89 


B0 


4.471 


4.50 


X 


X 




6 1 OriD 


6.700 


1.79 


B0.5 


4.420 


4.35 








NU Ori 


6.840 


2.09 


Bl 


4.383 


4.33 


X 


X 


X 


2 OriB 


6.410 


0.73 


Bl 


4.383 


3.96 








Par 1772 


8.370 


1.47 


B2 


4.294 


3.27 






X 


JW660 


9.620 


2.44 


B3 


4.272 


3.11 








9 1 OriB 


7.260 


0.49 


B3 


4.272 


3.28 




X X 


? 


'Weak- Wind' Sample 


6 2 Ori C 


8.240 


0.92 


B5 


4.140 


2.79 


? 




X 


6 l OriE 






B5 


4.140 


2.42 








6 1 OriF 






B8 


4.053 


1.84 








JW197 


10.16 


0.85 


B9 


4.025 


1.69 








JW153 


8.920 


0.13 


B9 


4.025 


1.89 






X 


JW831 


9.470 


0.00 


AO 


3.993 


1.54 








JW531 


10.55 


1.56 


Al 


3.975 


1.71 






X 


JW108 


10.19 


2.07 


A3 


3.940 


2.01 






? 


JW608 


11.89 


1.29 


A5 


3.917 


1.00 








JW165 


13.55 


4.20 


A7 


3.897 


1.49 








JW599 


19.20 


10.7 


A9 


3.867 


1.86 






? 



a Crosses indicate variable radial velocity or confirmed spectroscopic binary, dashes 
indicate constant radial velocity (Frost et al. 1926; Abt et al. 1991; Morrell & Levato 
1991). 

b Crosses indicate visual companions identified in speckle surveys (Weigelt et al. 
1999; Preibisch et al. 1999; Petr et al. 1998; Simon et al. 1999; Mason et al. 1998), 
dashes indicate negative results. 

c Optical binaries identified by eclipses are indicated by crosses (Antokhina et al. 
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1989). 

d Crosses mark HAeBe candidates (The et al. 1994; Tovmassian et al. 1997). '?' 
mark stars with L-band excess in the near-IR color magnitude diagram (Figure 11). 
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Table 2. COUP X-ray properties of early-type stars in the ONC 



COUP 


ID 


SpTy 


Offaxis 


A ix 


NetCts 


PSF 


Exp 


HR 


FBG02 


# 








a 




frac 


ks 




# 




'Strong- Wind' Sample 


809 


6 1 Ori C 


07 


0.32 


0.05 


63444 


0.01 


834.7 


26 ± 00 


542 


1232 


9 2 Ori A 


09.5 


1.95 


0.15 


3525 


0.01 


827.6 


0.09 ±0.01 


828 


745 


^OriA 


BO 


0.53 


0.13 


21283 


0.02 


834.7 


0.34 ±0.01 


498 


869 


6 1 OnB 


B0.5 


0.39 


0.12 


7901 


0.87 


832.9 


0.04 ±0.00 


584 


1468 


NU Ori 


Bl 


8.41 


0.40 


2237 


0.67 


728.6 


0.03 ± 0.00 


996 


1360 


# 2 OriB 


Bl 


2.69 


0.15 


1293 


0.86 


767.5 


0.21 ±0.02 


916 


349 


Par 1772 


B2 


4.58 


0.08 


127 


0.87 


758.6 


0.18 ±0.05 




1116 


JW660 


B3 


2.13 


0.07 


52865 


0.87 


818.7 


0.65 ±0.01 


746 


778 


9 1 Ori BE 


B3 


0.60 


0.13 


1506 


0.37 


832.9 


0.17 ±0.02 


519 








'Weak-Wind' Samp] 


e 








1473 


# 2 OriC 


B5 


3.92 


0.17 


26525 


0.86 


763.9 


0.38 ±0.01 


995 


732 


^OriE 


B5 


0.59 


0.11 


37523 


0.04 


834.7 


0.34 ±0.01 


495 




9 1 OnF 


B8 






< 29 




822.5 
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JW197 
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4.50 


0.18 


5311 


0.87 


811.7 


0.29 ±0.01 


103 


100 


JW153 


B9 


8.61 


0.38 


679 


0.89 


594.2 


0.14 ±0.02 


70 


1415 


JW831 


AO 


3.90 


0.08 


866 


0.87 


801.1 


0.08 ±0.01 


952 




JW531 


Al 






< 20 




815.9 








JW108 


A3.5 






< 66 




454.9 








JW608 


A5 






< 40 




801.2 






113 


JW165 


A7 


9.44 


0.32 


6724 


0.89 


212.2 


0.32 ±0.01 


81 


995 


JW599 


A9 


6.20 


0.16 


214 


0.87 


767.5 


0.75 ±0.13 


670 
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Table 3. COUP X-ray spectral parameters of early-type stars in the ONC 



COUP 


ID 


SpTy 


log7V H 


kT\ 


kT 2 


kT 3 


log EMi lo 


g EM 2 log EMz 


Xi (dof) 


log L t 


log Lt, c 


log (Lt,c/L bo i) 


# 






cm 


keV 


keV 


keV 


— 3 

cm ° 


zm ° cm ° 




erg s _1 


erg s _1 
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1.10 ( lo) 


90 c; 

zy.o 


on 

zy.o 


v n 

- f .y 


349 


Par 1772 


B2 


21.7 


0.44 


2.58 




51.9 


51.2 


1.04 ( 13) 


28.4 


28.9 


-7.9 


1116 


JW660 


B3 


21.7 


0.97 


3.33 




53.8 


54.3 


1.85 ( 263) 


31.2 


31.4 


-5.3 


778 


9 1 OriBE 


B3 


21.4 


0.63 


2.26 




53.0 


52.7 


0.70 ( 16) 


29.9 


30.1 


-6.7 
















'Weak-Wind' 


Sample 










1473 


6» 2 OriC 


B5 


21.4 


0.86 


2.86 




53.4 


53.9 


1.88 ( 166) 


30.9 


31.0 


-5.4 


732 


e 1 OriE 


B5 


21.1 


0.64 


1.58 


>15 


54.8 


55.3 54.2 


2.58 ( 190) 


32.3 


32.4 


-3.6 


142 


JW197 


B9 


21.0 


0.81 


3.06 




52.6 


53.0 


1.08 ( 61) 


30.1 


30.2 


-5.1 


100 


JW153 


B9 


20.8 


0.66 


1.39 




52.1 


52.2 


0.83 ( 29) 


29.3 


29.4 


-6.1 


1415 


JW831 


AO 


21.3 


0.26 


1.00 




52.6 


52.2 


1.11 ( 31) 


29.2 


29.6 


-5.6 


113 


JW165 


A7 


21.4 


0.87 


2.41 




53.3 


53.7 


1.54 ( 80) 


30.7 


30.9 


-4.2 


995 


JW599 


A9 


22.0 


1.89 






52.2 




0.71 ( 25) 


28.8 


29.2 


-6.3 
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Table 4. X-ray variability of strong-wind (SW) and weak-wind (WW) stars in the ONC 



COUP # 


ID 


SpTy 


Prob KS 


Prob MLB 


N F 


^char 
^total 


Ft i 

^char 






'Strong- Wind' Sample 








809 


6 1 OnC 


07 


V 


V 





1.0 


1.2 


1232 


2 OriA 


09.5 


V 


V 





1.0 


1.8 


745 


9 1 OriA 


BO 


V 


V 


2 


0.9 


2.9 


869 


6 1 OriD 


B0.5 









1.0 


1.0 


1468 


NU Ori 


Bl 









1.0 


1.0 


1360 


6 2 OriB 


Bl 


V 


V 


i 


0.9 


4.2 


349 


Par 1772 


B2 


V 


V 





1.0 


1.4 


1116 


JW660 


B3 


V 


V 





0.9 


1.5 


778 


6 1 Ori BE 


B3 




V 


1 


1.0 


3.3 






Weak- Wind' Sample 








1473 


9 2 Ori C 


B5 


V 


V 


3 


0.7 


3.0 


732 


6 1 OriE 


B5 


V 


V 


1 


0.9 


1.9 


142 


JW197 


B9 


V 


V 


4 


0.8 


51.3 


100 


JW153 


B9 


V 


V 


2 


0.9 


11.8 


1415 


JW831 


AO 


V 


V 





1.0 


1.6 


113 


JW165 


A7 


V 


V 


1 


0.9 


3.8 


995 


JW599 


A9 




V 





1.0 


1.6 
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Table 5. Comparison of flares on SW, WW and solar-type stars in the ONC 



Parameter 




Sample 




SW+WW WW 


1 M 




9 


7 


28 


^Vf, total 


15 


11 


41 


-^f,spec 


12 


9 


37 


Median duration (ks) 


13 


18 


65 


Median logLp (erg s -1 ) 


32.3 


31.2 


30.8 


Median logE^ (erg) 


36.2 


35.4 


35.5 



Note. - Numbers for 1 M stars from Wolk et al. 
(2005); N* = number of stars; -/Vf total = total number 
of flares; Nf tSpec = number of flares with spectral infor- 
mation 



Table 6. Periodic variability of SW and WW stars 



Name 


-Popt 


Px 


FAP 


Median 


Ampl 




d 


d 




cts ks -1 


cts ks -1 


9 1 OriC 


15.42 


15.51 


0.000 


72.6 


93.9 ±4.2 


JW660 


6.15 


6.67 


0.001 


61.3 


30.7 ±6.6 


JW165 


5.77 


4.61 


0.001 


7.8 


6.1 ±3.8 
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Table 7. Multiple systems with early-type primary in the COUP field of the ONC 



Name 


Type 


Comp 


Component separation 
Sep (mas) PA (°) 


Epoch 


Masses (M ) 
Pri Sec 


Ref. 


COUP 

# 


^OriC 


VB 


C1-C2 


38 ± 2 


208 


2001 


45 


> 6 


a 


809 


9 2 OriA 


VB 


A1-A2 


383 ± 10 


z91 


1997 


OK 

zo 


~ I 


b 


1232 




CD 
OJJ 


A 1 A "i 
Al-Ao 


100 






25 


~ 9 


b 




6 1 OriA 


VB 


A1-A2 


215 ±3 


357 


2001 


16 


~ 4 




745 


6 1 OriA 


eel 


A1-A3 








16 


~ 3 


c 




NU Ori 


VB 


A-B 


471 ± 17 


98 


1997 


14 


~ 1 


b 


1468 


NU Ori 


SB 


A-C 


80 






14 


~ 3 


b 




9 1 OriB 


VB 


B1-B2 


942 ± 10 


255 


2001 


7 


~ 4; 1.6 


a b 


778+766 


9 1 OriB 


VB 


B2-B3 


117 ±4 


210 


2001 




~ 3; 0.7 


a b 




6 1 OriB 


VB 


B1-B4 


609 ±8 


298 


2001 


7 


~ 0.2 


a b 




1 OriB 


SB 


B1-B5 


30 






7 


~ 2.7 


bd 




9 2 OriC 


SB 


C1-C2 












e 


1473 



Note. — VB = visual binary, SB = spectroscopic binary, Comp = component, Pri = primary, 
Sec = secondary, PA = position angle (degrees East of North) 

a Schertl et al. (2003) 

b Preibisch et al. (1999) 

c Lloyd & Strickland (1999) 

d Palla & Stahler (2001) 

e Corporon & Lagrange (1999) 



-42- 



Table 8. Tentative classification of early- type ONC X-ray sources 



COUP # 


Name 


SpTy 


Luminosity Spec 


Var 


Disk 


Magn 








'Stron 


5- Wind' Sample 








809 


^OriC 


07 


>W 


S 




RM,B 


1232 


# 2 OriA 


09.5 


W W 


w a 






745 


^OriA 


B0 


>W 


T 




NT 


869 


^OriD 


B0.5 


W W 


W 






1468 


NU Ori 


Bl 


W W 


w 


flR 




1360 


<9 2 0riB 


Bl 


w 


T 






349 


Par 1772 


B2 


w 




flR 


CP? 


1116 


JW 660 


B3 


>w 


s 




RM 


778 


9 1 On BE 


B3 


w 




nIR 




'Weak-Wind' Sample 


766 


6 1 On BW 


? 




T 






1473 


e 2 Ori c 


B5 


T 


T 


flR 




732 


^OriE 


B5 


T 


T 




NT 




^OriF 


B8 










142 


JW 197 


B9 


T 


T 






100 


JW 153 


B9 


T 


T 


flR 




1415 


JW 831 


AO 


T 










JW 531 


Al 






flR 






JW 108 


A3 






nIR 






JW 608 


A5 










113 


JW 165 


A7 


T 


T/S 




RM 


995 


JW 599 


A9 


T 




nIR 





Note. - Classes: W = hydrodynamic wind (many small shocks), S 
= smoothly variable lightcurve, T = T Tauri type emission, nIR = near- 
infrared excess, flR = far-infrared excess, RM= rotational modulation, B = 
detected magnetic field, CP = chemically peculiar, NT = non-thermal cm 
emission. 

a This star exhibited powerful flares during the earlier Chandra observa- 
tion described by Feigelson et al. (2002), and thus might also be classified 



"T" here. 



43 - 



-44- 




Log T,„ [K] 

Fig. 1. — Hertzsprung- Russell diagram with pre-main sequence models from Siess et al. 
(2000). Dotted grey curves show isomass evolutionary tracks, solid grey curves show 
isochrones, and solid black jagged contours indicate the size of the convective envelopes 
labeled in units of the stellar radius. Symbols show spectroscopically confirmed ONC stars 
observed in the COUP: triangles are magnetically active (MA) stars, and circles are weak 
wind (WW) stars. The strong wind (SW) stars are located outside the plotted range to the 
left. Grey symbols indicate stars not detected in the COUP. 
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Fig. 2. — Relation between the emission measures from the two components of the spectral 
fits from Table 3. For 3-temperature models, we plot the EM corresponding to the two 
lower temperatures which dominate the emission. The region included within the dotted 
lines marks differences of a factor two and lower. Squares are SW stars, and circles are WW 
stars. Three outliers are labelled by their COUP source numbers. 
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Fig. 3. — Ratio of the maximum to minimum count rate blocks in the MLB test for all sources 
in the COUP optical sample plotted versus the effective temperature. If the maximum count 
rate corresponds to a flare the plotting symbol is filled. The thick grey lines mark the median 
values for these ratios in each class: MA = magnetically active, WW = weak-wind, and SW 
= strong-wind. 
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Fig. 4. — Fraction of COUP sources which are X-ray variable as a function of source bright- 
ness. The curves show the distribution for MA stars in the COUP according to Kolmogorov- 
Smirnov and Maximum Likelihood Blocks criteria. The individual hot SW and WW stars 
are plotted as open boxes and circles, respectively. 
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Fig. 5. — Lightcurves of COUP count rate and hardness ratio for detected SW and WW 
stars in order of decreasing spectral type. Segments derived from the MLB analysis are 
overplotted with horizontal lines onto the binned lightcurve. The hardness ratios have been 
evaluated in individual segments. 
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Fig. 5. — continued 
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Fig. 6. — Histogram showing the distribution of flare durations for the 11 flares on WW 
stars. The exponential 0.1 x e ( ~*/ 100fcs) from Wolk et al. (2005) is overplotted. 
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Fig. 7. — Folded X-ray lightcurves and hard- 
ness ratios for the three stars with optical 
photometric periodicities. 6 1 Ori C is folded 
with the optical period using the ephemeris 
by Stahl et al. (1996). JW660 and JW165 
are folded using the period derived from the 
X-ray data. Data from consecutive cycles are 
shown with different plotting symbols. The 
hardness ratio is examined in 5 phase bins. 
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Fig. 8. — COUP image of the region near the 6 1 Ori B quintet displayed at very high res- 
olution (0.125" pixel size). The yellow and magenta x symbols show the positions of the 
optical components. 
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Fig. 9.— ACIS spectra of 9 1 Ori B East (SpT B3; COUP 778) and 9 1 Ori B West (?; 
COUP 766). The data are shown with error bars, the continuous histogram is the best- 
fit model from XSPEC, and the bottom panels show residuals between the data and models. 
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Fig. 10. — COUP image showing the X-ray detection of the B1.5 star Par 1772 (red cir- 
cle), resolved for the first time in X-rays from adjacent brighter X-ray sources. 'C-XXX' 
denotes the COUP source number, and the green circle shows the ROSAT HRC error circle 
previously associated to Par 1772. 
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Fig. 11. — Near-IR color-color diagrams for SW and WW stars in the COUP field. Smaller 
symbols are for stars not detected in the COUP, and filled symbols for known or suspected 
binaries. Plotting symbols indicate the following properties: triangles are JiM.S'-excess stars 
from The et al. (1994), squares are chemically peculiar (CP) star candidates from Renson 
et al. (1991), and circles are the remaining normal stars. The solid and the dotted curves are 
the loci of dwarf and giant stars from Bessel & Brett (1988), and the dashed lines indicate 
a reddening corresponding to Ay = 10 mag using the extinction law by Rieke & Lebofsky 
(1985). 
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Fig. 12. — X-ray efficiency as a function of effective temperature. The dotted lines indicate 
the threshold values defined in Sect. 7.1. The plotting symbols are defined in Figure 11. 



